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Fully microscopic RPA/LDA calculations of the dipole 
plasmon for very large neutral and charged sodium clusters, 
Na^ + , in the size range 950 < N < 12050 are presented for 
the first time. 60 different sizes are considered altogether, 
which allows for an in-depth investigation of the asymptotic 
behavior of both the width and the position of the plasmon. 
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The surface dipole plasmon in large metal flusters is 
currently attracting .significant experimentaltTH and the- 
oretical attentionpti e.g., a time-resolved analysisu of 
second harmonic generation via femtosecond pump-probe 
studies was recently successful in establishing an upf 
limit of about 25 nm for the validity of the l/i?-lawB 
(where R is the cluster radius) associated with the width 
of the dipole plasmon in large Na^v clusters. On the 
theoretical side, this Landau-damping-type width, T = 
AUvf/R, in large metal clusters has been-derived repeat- 
edly with various analytical approachesQtl (which nec- 
essarily utilize simplified approximations for the dipole 
matrix elements of the residual Coulomb force), but a 
microscopic justification of it is still lacking. In gen- 
eral, detailed microscopic studies (the most systematic 
among them using the matrix random phase approxi- 
mationilpcal density approximation (matrix- RPA/LDA) 
versiontrli3 of the linear response theory) have been resc- 
tricted-tp rather small clusters in the size range 8 < N < 

338.yy 

In the absence of systematic matrix- RPA/LDA cal- 
culations above N = 338, various other theoretical ap- 
proaches have been most recently introduced with the 
expressed purpose of providing a framework for micro- 
scopic studies of the dipole plasmon in large metal clus- 
ters; they include the in-real-time version of the time- 
dependent local density approximation the modified- 
with-separable-residual- forces RPA-, j and the separation- 
of-collective-coordinates approach£l Due to the compu- 
tational effort involved, however, these methods were in 
their turn restricted to a few (about 10) smaller sizes, 
in the rangeO of 8 < N < 400 (or even in the rangeD 
8 < N < 950), and as a result they were unable to in- 
vestigate the asymptotic behavior of the position and the 
profile of the plasmon. In addition, through an identifica- 
tion of the width with the variance of the RPA response, 
the authors of the recent Rapid Communications came 
via a semianalytical analysis to the conclusion that the 
plasmon width should exhibit an 1/vi? asymptotic de- 



pendence, in sharp contrast with the l/R law. 

In this paper, with the help of the new generation of 
powerful computers available today, I present systematic 
matrix- RPA/LDA investigations of the dipole plasma ex- 
citation for neutral and charged sodium clusters in the 
size range 950 < N < 12050. In each_case 30 different 
sizes were calculated]!.] which allowedlij for the deter- 
mination of the asymptotic size dependence of both the 
position and the profile of the plasma resonance. 

The theoretical aspects of the matrix- RPA/LDA ap- 
proach have been presented in Ref. [Ic]. Here I only dis- 
cuss briefly how the matrix-RPA equations (see Eq. (5) 
of Ref. |ToD can be reduced to an equivalent problem with 
half the dimensions of the original matrix (such a reduc- 
tion of dimensions is instrumental for the successful im- 
plementation of the matrix-RPA/LDA algorithm in the 
case of the very large metal clusters considered here). 
For such a reduction, it is sufficient to add and subtract 
the top and bottom rows of the original-matrix equation. 
One then finds the equivalent equationjlj 
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where the matrix A is given by Eq. (14) of Ref. A 
in general is the multipolarity of the excitation (A = 1 
for the case of the dipole), £„ are the RPA eigenvalues, 
and D is a diagonal matrix depending exclusively on the 
unperturbed single particle-hole energies, i.e., D p i ltP >h' — 
{e p — eh)5 pp '5hh' ■ The usual forward- and backward-going 
RPA amplitudes W" = X u and W v _ = Y v , which are 
needed in order to calculate the oscillator strengths f v 
(see section II. C of Ref. [!(]), are given through the eigen- 
vectors Z v as follows: 



(2) 



where the proportionality constant C is determined via 
the RPA-cigenvector normalization condition. 

For each cluster Na^ + , the microscopic RPA calcula- 
tion yields a discrete set of oscillator strengths /„ (nor- 
malized to unity) associated with the RPA eigenvalues E v 
(see inset of Fig. 1). By folding the oscillator strengths 
with Breit-Wigner shapes normalized to unity, one can 
calculate a smooth photoabsorption cross section, a, per 
valence electron as follows: 

a(E) = 1.0975 (eV A 2 ) V f v P B w(E; E v , lv ) , (3) 



where 7„ denotes the intrinsic widths of the auxiliary 
Breit-Wigner profiles, Pbw- 
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Fig. 1 displays the calculated a's for Nag52 and Nai 2 o68, 
namely for the smallest and the largest size considered in 
the present paper (intrinsic widths of 7^ = 7 = 0.16 and 
0.076 eV were used for Nag52 and Nai2068, respectively). 
It is apparent that the FWHM (full width at half max- 
imum) is smaller for Nai2068 (~ 0.147 eV) compared to 
the case of Na952 (~ 0.319 eV) and that the maximum 
of the photoabsorption cross section for N — 12068 is 
blueshifted with respect to the case of N = 952. 

In addition to the photoabsorption profiles, Fig. 1 also 
displays the variance £ (horizontal solid bar) associated 
with the RPA dipole response in the .case of Nai2068 
(in the RPA the variance is calculated^ via the rela- 
tion E 2 = (111,2/1710) — (mi/mo) 2 , where the moments nik 
are given by rrik = £^~^f v ). Even with the overesti- 
mation of the FWHM due to the folding procedure (see 
below), it is seen that the variance is substantially larger 
than the corresponding FWHM in the case of Nai2068! 
thus the variance should not be used as a substitute fps. 
the actual width of the dipole plasmon, as was doneH 
in Ref. |^. Furthermore, to correct for the overestimation 
due to the folding, one needs to subtract the intrinsic 
width 7„ from the original FWHM (see Fig. 10 of Ref. 
|l0| and accompanying discussion) . With this correction, 
the plasmon widths are ~ 0.16 eV and ~ 0.07 eV for 
Nag52 and Nai2068: respectively; these values are even 
smaller than the uncorrected FWHM's compared to the 
corresponding variances. 

The method of folding described above has been used 
in most of the earlier publications in order to extract the 
FWHM of the photoabsorption profiles. This method, 
however, becomes cumbersome in the present study 
where the size changes substantially going from the small- 
est to the largest cluster. Indeed, smaller 7^'s must be 
used the larger the cluster size, and such a procedure is 
difficult to be determined uniquely. Instead of the fold- 
ing method, a more uniform procedure is needed. Using 
the fact, that for large cluster sizes it can be analytically 
shownO that the RPA resp&iase itself tends asymptoti- 
cally to a well concentratedE3 distribution of oscillator 
strengths /„ exhibiting a Breit-Wigner profile, i.e., 

f v ^v/[{£ v -£f + {T/2f], (4) 

this paper has opted for a method based directly on the 
distribution of the /i/s; namely, the Landau-damping- 
type width r is determined by the minimum interval 
around the RPA energy centroid £ which contains 50% of 
the total oscillator strength, i.e., V is such that J2i fi = 
0.5 with £i e (£ - T/2,£ + T/2). The value of 50% 
used here is precisely the percent of the total oscillator 
strength contained within the energy interval associated 
with the FWHM of the single Breit-Wigner profile [see 
Eq. (Q)], and the RPA centroid is given by £ = mi/m . 

Fig. 2(a) displays the calculated (according to the lat- 
ter method) RPA widths (r, open circles) and variances 
(£, open squares) for Na^ clusters in the size range 
952 < N < 12068 as a function of the inverse radius 



(1/i?) of the jellium background (R — rgN 1 / 3 , and the 
Wigner-Seitz radius for sodium was taken as r s — 4 a.u.). 
Fig. 2(b) displays the same quantities as a function of 
1/y/R. Before proceeding further, it needs to be no- 
ticed that the corrected FWHM's for Nag52 and Nai2068 
(namely the values 0.16 and 0.07 eV) according to the 
folding method are in good agreement with the widths T 
derived from the 50% method. 

A first observation is that the variances are systemat- 
ically larger than the widths (see also Fig. 1). Second, 
it can be seen that the widths are not proportional to 
1/VR, since the associated straight line does not pass 
through the origin of the axes [Fig. 2(b)]. On the con- 
trary, in the case of Fig. 2(a), the straight dashed line 
that passes on the average through the points denoting 
the widths definitely passes through the origin of axes. 
This provides an unequivocal proof that T cx 1 / R. No- 
tice that for the smaller sizes the precise values of the 
widths exhibit noticeable scattering around the average 
line specified by the 1/i? law; however, such scattering 
becomes progressively smaller the larger the cluster sizes. 
The solid line in Fig. 2(a) denotes the T — Afivp/R 
law according to the analytical result of Ref. ||, namely 
whenEil A = A ana = 0.46. For the numerical RPA calcu- 
lation, one finds A num = 0.58, which is close to A ana . The 
difference between the two slopes is due to the infinite- 
well approximation (and ensuing neglect of electronic 
spillout) which was invoked during the analytical deriva- 
tion. To further test the 1/R dependence, I have also 
plotted in Fig. 2(a) the uncorrected FWHM's of the pho- 
toabsorption profiles for Nag52 and Nai2068 (solid dots). 
Again, and in sharp contrast with the variances, even 
these uncorrected FWHM's follow clearly an 1/R depen- 
dence. 

Finally, _it_ has been derived through semiclassical 
argumentaj'tH that the size dependence of the variance 
should be propotional to X/yR. Fig. 2(b) does not indi- 
cate such a relation; much larger sizes are needed before 
the microscopic results for E do converge to the expected 
1 / \R relation. 

Fig. 3 portrays the size evolution of the quantities as- 
sociated with the position of the RPA surface plasmon, 
namely, it displays the quantities £\ = (mi/m-i) 1 / 2 , 
£3 = (ma/mi) 1 ' 2 , and the centroid £ as a function of 
the inverse cluster radius. Fig. 3(a) corresponds to the 
case of neutral Najv clusters, while Fig. 3(b) displays 
the case of multiply cationic Na]^ + clusters. In all in- 
stances, the calculated RPA values honor the theoreti- 
cally expectedEj inequality £1 < £ < £3, but the cen- 
troid lies much closer to £\ than to £3. For the neutral 
clusters, it is seen that the calculated points for £\ and 
£ lie closely on two straight lines, which cross the energy 
axis (R = 00) at about 3.403 eV, in very good agreement 
with the expected classical Mie limit for sodium clusters, 
i.e, huJMic = ^buik/v^ = 3.4 eV (for r s = 4 a.u.). From 
the relation, £ — 3.403 eV — rj/R, one can determine 
77 = 4.68 (eV a.u.). The calculated points for £3 do not 
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fall on a straight line crossing the energy axis at about 
3.4 eV, although such a behavior hag, been predicted from 
semiclassical sum-rules arguments jl3'B2l apparently much 
larger sizes are needed before the size dependence of £3 
converges to this expected asymptotic behavior. 

Comparing the RPA results in Fig. 3(b) with corre- 
sponding results in Fig. 3(a), one sees that the position 
of the dipole plasmon in the case of charged clusters is 
blueshifted with respect to that of the neutral clusters (all 
three quantities £\, £, and £3 exhibit this trend). The 
extent of this blueshift depends on the ratio Z/N, and 
naturally for the larger charged clusters the calculated 
positions tend to converge to the corresponding results 
for the neutral clusters; for Z = 10, this convergence is 
definitely recognizable for the largest sizes studied here. 

Finally, Fig. 4 displays the size evolution of the RPA 
widths r and variances E of the charged Na]^ + clusters 
as a function of 1/R. As was the case with the position, 
the charge state influences the values of T and E; in par- 
ticular a reduction in magnitude of these quantities can 
be seen compared to the case of neutral clusters. This 
reduction depends again on the ratio Z/bL and it is more 
pronounced the smaller the cluster size.E-3 

In conclusion, fully microscopic RPA/LDA calcula- 
tions for neutral Najv and charged Na]^ + (950 < N < 
12050) clusters were presented. In each case, some 30 
different sizes were considered, which allowed for an in- 
depth investigation of the asymptotic behavior of both 
the position and the width of the dipole plasmon. In par- 
ticular, it was found that asymptotically the RPA width 
becomes proportional to the inverse cluster radius, and 
that for neutral clusters this trend is already well devel- 
oped within the size range considered here. 
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formed at the GIT Center for Computational Materials 
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No such closed semiclassical asymptotic behavior has been 
derived for Si and S. 
24 A similar dependence of the plasmon width (and of the 
position as discussed earlier) on the charge state was also 
found in the case of small clusters with N — Z < 40 (see 
Refe.|0@). 

FIG. 1. RPA photoabsortion cross sections for Nag52 
(dashed line) and Nai2068 (solid line). The solid horizontal 
bar denotes the RPA variance E in the case of Nai2068- Inset: 
the oscillator-strength distribution for Nai206S- 

FIG. 2. RPA widths P (open circles) and variances E 
(open squares) of the dipole plasmon for neutral Najv clus- 
ters (952 < N < 12068), as a function of (a) the inverse 
cluster radius, 1/R and (b) the inverse of the square root of 
the radius, l/y/R. The solid line in (a) corresponds to the 
analytical result for F according to Ref. ^. The solid dots 
in (a) denote the uncorrected FWHM of the photoabsorption 
profiles (see Fig. 1) for Nag52 and Nai2068 generated with the 
folding method. 

FIG. 3. The RPA positional quantities Si (open squares), 
S (centroid, solid circles), and S3 (open triangles) of the dipole 
plasmon, as a function of 1/R for (a) neutral Najv and (b) 
charged NaJ^ + clusters. The straight line in (b) is a duplicate 
of the lowest line in (a) and was drawn as a guide to the eye. 



FIG. 4. RPA widths F (open circles) and variances E 
(open squares) of the dipole plasmon for charged NaJ^ + clus- 
ters (962 < N < 12078), as a function of the inverse cluster 
radius, 1/R. The straight lines are duplicates of dashed lines 
in Fig. 2(a) and were drawn as guides to the eye. 
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